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We present ultrafast x-ray diffraction experiments on thin films of metallic SrRuC>3 (SRO) after 
their excitation with ultrashort intense laser pulses. Depending on the layer thickness, the data 
exhibit a transient splitting of the (002) SRO Bragg peak evidencing the generation and propagation 
of sharp acoustic strain waves. These distinct structural dynamics are due to the exceptionally fast 
electron-phonon relaxation that gives rise to a quasi-instantaneous thermal stress in SRO. The 
interpretation is corroborated by numerical simulations which show excellent agreement with the 
experimental findings. Despite the qualitatively different lattice dynamics for different SRO layer 
thicknesses, we identify a universal evolution of the transient average layer strain. The inferred 
discrepancy of the thermal stress profile from the excitation profile may hint toward a temperature- 
dependent effective Griineisen parameter of SRO. 

Keywords: Ultrafast, X-ray diffraction, Thin films, Lattice dynamics 



I. INTRODUCTION 

In general, the internal energy of condensed matter is 
spread over various degrees of freedom such as electrons, 
lattice, magnetization and polarization. Intense and ul- 
trashort laser pulses can directly excite one or more of 
these degrees of freedom within a given material. The 
subsequent dynamics of the individual subsystems and 
thus the entire system are then governed by the coupling 
strengths of the different degrees of freedom. For in- 
stance, structural dynamics of a crystalline lattice can be 
induced directly by infrared/THz absorption 1 - and Bril- 
louin/Raman scattering 2 -! 3 - or indirectly by an initial elec- 
tronic excitation (interband or intraband 4 ^ 5 -) and subse- 
quent electron-phonon coupling 5,6 . Depending on the 
material of interest, the electronic and phononic sub- 
systems may also be coupled to other degrees of free- 
dom such as polarization 6 ' 7 , magnetization 8 - - — or charge 
and/or orbital order— ~—. The understanding of these 
complex physical processes for a given material is of 
fundamental interest and, moreover, promises techno- 
logical advances in the fields of signal processing, data 
storage and sensors as well as novel x-ray optics for ul- 
trafast studies 1 ^—. In particular, the coupling of vari- 
ous degrees of freedom is often mediated by the lattice. 
This triggered an enormous interest in time-resolved scat- 
tering techniques during the last two decades in order 
to monitor the laser-induced changes of the structural 
propertie s 4 ' 18- — . 

A material of particular interest is the "bad metal" 
SrRu03 (SRO) due to its various complex physical prop- 
erties such as itinerant ferromagnetis m 22 : 23 , ultrafast 
magnetostriction and electron-phonon coupling 9 -^, neg- 



ative spin polarization^ 5 -, orbital ordering^ 6 - and non- 
Fermi liquid behaviou r 2 i 28 . This material can be epitax- 
ially grown on single-crystal substrates (such as dielectric 
SrTi03 [STO] ) with high structural perfection 2 ^. In com- 
bination with a fast electron-phonon relaxation of 200 fs 
or less 9 -^ and a very high damage threshold 30 , SRO is 
perfectly suited as transducer material for the generation 
of coherent longitudinal acoustic (LA) phonons 3 - 1 -. 

This report focuses on the structural dynamics of lay- 
ered crystalline solids induced by ultrashort laser pulses. 
In particular, we discuss the generation and evolution of 
acoustic deformations of SRO thin films on a supporting 
STO substrate. We utilize the experimental method of 
ultrafast x-ray diffraction (UXRD). This technique em- 
ploys the pump-probe scheme in which the laser-induced 
structural dynamics are probed by an ultrashort hard x- 
ray pulse at different time delays r after the arrival of 
the excitation (pump) pulse. After a brief introduction 
of the theoretical framework which describes the ultra- 
fast build-up of laser-induced thermal stress in SRO, we 
present results of UXRD experiments on two thin films 
with thicknesses smaller and larger than the optical pen- 
etration depth of the 800 nm pump light. The tran- 
sient changes in the UXRD data readily evidence a com- 
plex formation and propagation of LA phonon wavepack- 
ets. However, qualitatively different features appear for 
the two different samples. In particular, the thicker 
SRO layer exhibits a splitting of the Bragg peak as op- 
posed to a continuous shift in case of the thinner film. 
The experimental data can be precisely simulated by 
means of numerical models for the photoinduced struc- 
tural dynamics 3 ^ and the dynamical diffraction of x-rays 
from these transient crystal structures 3 - -. The universal 



2 



features of the lattice dynamics are analyzed in detail by 
considering the spatiotemporal strain fields and the po- 
tential and kinetic energy of the photoexcited thin film 
as well as the substrate. 



II. HEATING OF SRO THIN FILMS BY 
ULTRASHORT LASER PULSES 

The topic of laser-induced heating and transport prop- 
erties of elemental metals (or metal layers) is a fairly well 
understood process^ 2 - and has been discussed many times 
in literature^—. In the following we want to give a brief 
summary and apply the standard theoretical models to 
the case of thin films of SRO on a STO substrate. 

The standard model for describing the transient pro- 
cesses in laser-heated metals is the two-temperature 
model (TTM) proposed by Anisimov et al?—. It assumes 
that the optical energy of the laser pulse is entirely ab- 
sorbed by the conduction band electrons in the metal. 
The electronic system then rapidly thcrmalizes towards 
an elevated electron temperature via electron-electron 
scattering processes. The temperature of the electronic 
system T e then differs from the temperature of the lattice 
T; (phononic system) and subsequent elect ron-phonon 
collisions transfer energy from electrons to phonons until 
the two subsystems reach thermal equilibrium. Due to 
lattice anharmonicities the incoherently excited phonons 
produce thermal stress which eventually leads to the ther- 
mal expansion of the metal. For most metals the elec- 
tron relaxation time r e is on the order of a few tens of 
femtoseconds (fs)2S and is thus shorter than the typical 
timescales of laser pulse durations and all other dynam- 
ical processes involved. In particular, this assumption 
holds for SRO which exhibits a very short electron re- 
laxation time of r e w 4.2 fs at 145 This validates 
the consideration of an electron temperature at all times 
and restricts the electronic heat transport to be diffu- 
sive. Typically, the linear dimensions of the excitation 
and probe area on the sample surface is much larger than 
the penetration depth or the film thickness of the metal. 
Under these circumstances the differential equations of 
the TTM can be restricted to one spatial dimension and 
read as follows: 

(i) 

where C e /i and k e /i are the electronic/lattice heat ca- 
pacity and conductivity, respectively, G is the electron- 
phonon coupling factor and S(x, t) is the heat source de- 
termined by the laser pump puls o 34 ' 35 . The usual consid- 
erations of the TTM for elemental metals epmloy the fact 
that the heat is dominantly conducted by the conduction 
band electrons and that the phononic heat conductivity 



is comparably small. This allows the omission of the 
first term in (J2]>22r— . In SRO, however, the heat is car- 
ried by the electrons and lattice in approximately equal 
part s 38 ' 39 which is why we keep this term in @. 

As motivated above, the electron-electron scattering 
rate is very large in SRO^i. The electrons are thus not 
able to ballistically transport energy out of the excited 
region into deeper parts of the sample. In addition, 
SRO is known to have a very fast electron-phonon re- 
laxation time (i.e. large G) on the order of a few hun- 
dred femtosecond*^ which is also much faster than any 
diffusion processes of electrons and phonons. When con- 
sidering the structural dynamics in thin SRO films on 
the timescale of a few tens of picoseconds, we can thus 
disregard the diffusion terms in ([T]) and {2J . This simpli- 
fies the problem considerably and the eventual expansion 
profile (caused by the thermal stress profile) in SRO is 
proportional to the exponentially decaying profile of the 
deposited energy density where the proportionality con- 
stant is given by the Griineisen parameter 7^. The en- 
ergy density profile is given by the derivative of Lambert- 
Beer's law and thus defined by the optical penetration 
depth £ opt of the pump laser light at 800 nm wavelength. 

An implication of the very fast electron-phonon cou- 
pling in SRO is the fact that the thermal stress is built 
up quasi-instantaneously which launches coherent LA 
phonons up to very high frequencie s 9 ' 17 ! 24 ' 30 " — . Similar 
to the temperature considerations above the structural 
dynamics are reduced to one spatial dimension. The 
coherent longitudinal (plane) strain waves traverse the 
excited metal layer until the entire coherent vibrational 
energy has left into the substrate and a quasi-statically 
thermally expanded layer is left. In the following sec- 
tions we show that the coherent phonon dynamics in- 
side the excited SRO thin film and the underlying STO 
substrate can be accurately monitored by UXRD exper- 
iments. For the later analysis we define the character- 
istic timescale of sound propagation through a thin film 
by T soun( j = d/uscmnd where d is the film thickness and 
W S ound is the longitudinal sound velocity perpendicular to 
the sample surface. 

III. EXPERIMENTAL RESULTS 

The UXRD experiments presented and discussed in the 
following were conducted on thin SRO films of different 
thickness epitaxially grown on a STO substrate. The 
samples were excited by ultrashort optical laser pulses 
at 800 nm wavelength and the triggered structural dy- 
namics within the first few tens of picoseconds were ob- 
served by UXRD employing the Plasma X-Ray Source 
(PXS) at the University of Potsdam, Germany 41 . This 
laser-based x-ray source utilizes hard x-ray pulses at the 
characteristic Cu K Q i and K Q 2 line (-Epxg = 8.048 keV 
and E ( pl s = 8.028 keV) with a pulse duration of w 200 fs. 
The PXS thus provides the temporal and spatial resolu- 
tion required for studying coherent acoustic phonon dy- 
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FIG. 1. Static XRD curves (symbols) and dynamical XRD simulation of a (a) 15.4 nm and a (b) 94.8 nm SRO layer on STO 
around the (002) reflections. The data were recorded at the ESRF (E ES rf = 12.0 keV) and at BESSY II (£ B essy = 8.66 keV), 
respectively, at x-ray photon energies differing from the characteristic Cu K Q i line (-Bpxs = 8.048 keV) generated by the PXS. 
The Bragg angle axes in panel (a) and (b) were thus converted to an artificial x-ray photon energy -Epxs to be comparable 
to the plots (c) and (d). The lower panels present the UXRD data recorded with the PXS on the SRO Bragg peak of the (c) 
thinner and (d) thicker SRO thin film. The excitation fluences of the 800 nm pump pulses are 30 mj/cm 2 and 20 mj/cm 2 , 
respectively. 



namics in thin crystalline films with a thickness below a 
few hundred nanometer. 

The SRO thin film samples were pre-characterized by 
static x-ray diffraction (XRD) at synchrotron-based x- 
ray sources. The results of the static 6-20 scans around 
the respective (002) reflections are represented by the 
symbols in Fig. [Ha) an d (b). As expected, the thicker 
SRO film exhibits a narrow and intense Bragg reflection 
whereas the Bragg peak of the thinner SRO film is much 
broader and weaker. The solid lines are results of dy- 
namical XRD simulations used for the determination of 
the layer thickness and c-axis lattice parameter. The 
simulations revealed a thickness d± — 15.4 nm and an 
out-of-plane lattice parameter of c\ = 3.9525 A for the 
thinner SRO layer (Fig. [Ha)) an d di — 94.8 nm and 
c 2 = 3.9493 A (Fig. [TJb)) for the thicker layer. The 
different lattice parameters are consistent if one consid- 
ers the epitaxy with the cubic STO substrate (c su b = 
3.905 A) and the relaxation of the substrate-induced 
tetragonally distorted SRO unit cell as the layer thickness 
increases. Employing the longitudinal sound velocity of 



SRO, w S ound = 6.312 nm/ps (Ref. 1421). the derived layer 
thicknesses imply sound transit times of T s ^ nd = 2.45 ps 

(2) 

and Tg 0l ; nd = 15.0 ps for the thinner and thicker SRO 
layer, respectively. The two samples are chosen in order 
to represent the limiting cases of layers with thickness 
smaller and larger than the literature value of the optical 
penetration depth at 800 nm, £ p t « 52 nm (Ref. [27]), 
respectively. 

The UXRD data recorded at the PXS on the thinner 
and thicker SRO layer using a pump fluence of 30 and 20 
mJ/cm 2 are shown in Fig.QJc) and (d), respectively. Due 
to the limited signal-to-noise ratio and angle-resolution of 
the PXS the UXRD data quality is poorer for the thinner 
SRO layer and the respective SRO Bragg peak of the un- 
excited sample is not clearly separated from the substrate 
peak. Nevertheless, the data can be unambiguously an- 
alyzed in terms of transient SRO Bragg peak shifts by 
proper substrate subtraction and thus allow for a clear 
understanding of the coherent and incoherent phonon dy- 
namics in the photoexcited SRO layers. 

The UXRD data exhibit transient changes of the SRO 
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FIG. 2. Normalized transient shift and splitting of the SRO 
Bragg peak extracted from the UXRD data shown in Fig.QJc) 
and (d). The continuous shift observed in the thinner SRO 
layer is represented by the blue squares (single- Gaussian fit- 
ting). The empty and solid black bullets show the transient 
Bragg angles of the decaying initial and the rising shifted 
peak for the thicker SRO layer (double- Gaussian fitting). For 
7"/Tsound > 1 the red diamonds indicate the continuous back- 
shift of the thicker SRO layer (single- Gaussian fitting). The 
intensity-weighted shift for the thin SRO layer is shown by 
the red half-filled bullets. The error bars represent the 68% 
confidence interval of the fitting parameter and lines are a 
guide to the eye. 



Bragg peak positions just after the laser excitation at 
t = 0. This directly implies that the quasi- instantaneous 
heating of the SRO lattice triggers certain structural dy- 
namics inside the SRO layers. At first sight, these dy- 
namics appear to be qualitatively different. In case of 
the thinner SRO layer we observe a continuous shift of 
the Bragg peak towards lower angles for < t < T soun d 
(region I), followed by a slight continuous backshift for 
^sound < t < 2T soun d (region II) until it reaches a new 
quasi-stationary position for r > 2T sounc i (region III). 
This quasi-stationary expansion represents the thermal 
expansion of SRO due to the absorbed energy of the ex- 
citing laser pulse. In contrast, the thicker SRO layer does 
not exhibit a continuous shift for < r < T soun( j (region 
I). Instead, we find a splitting of the initial Bragg peak 
into two distinct reflections at intermediate times until 
the initial peak has disappeared. Similar Bragg peak 
splittings have been previously observed in photoexcited 
bulk crystals and thin films^ i 43 ' 44 , however, either the 
signatures were relatively weak compared to the bulk re- 
flection or a thorough description of the underlying struc- 
tural dynamics is missing. The comparison of the exper- 
imental results in Fig. [TJc) and (d) rises the question 
whether the structural dynamics responsible for the ob- 
served features are indeed qualitatively different. This 
issue is addressed in the subsequent sections. 

We extracted the transient Bragg angles of the mea- 



sured SRO peaks by fitting the data with single and 
double- Gaussian functions. The obtained peak positions 
are displayed in Fig. [2] where we plot the Bragg angle 
change normalized to the quasi-stationary value after 
2T soun d, AO(T)/AO(2T soun d), versus time delay in units 
of Tsound- Since the angle changes are fairly small they 
are proportional to the average strain variations of the 
SRO layer. This is indicated by the secondary y-axis in 
Fig. [5] presenting the normalized total SRO strain. The 
plot verifies the features visible in the contour plots in 
Fig- [He) and (d). The thinner sample exhibits a contin- 
uous shift of a single SRO peak (blue squares). In con- 
trast, the SRO peak in the thicker sample for normalized 
time delays up to T soun d (region I) shows a peculiar split- 
ting. The Bragg angles of the initial and displaced peaks 
are represented by the black empty and solid bullets in 
Fig. [21 respectively. In addition, we observe that both in- 
dividual peaks shift to higher Bragg angles as they decay 
and rise. Despite the unequal behaviour of the respective 
SRO peaks in region I, the normalized transient shifts in 
region II and III (r > T soun d) appear to be comparable in 
both SRO layers. In addition, we evaluate the intensity- 
weighted transient shift of the thicker SRO layer in region 
I which is shown by the red half-filled bullets in Fig. [5] 
This curve represents the sum of the individual shifts 
weighted by the intensity of the respective Bragg peak at 
each time step. We find that this curve coincides with 
the transient peak shift of the thinner SRO layer. This 
immediately implies that the time-evolution of the aver- 
age strain is identical in both layers as will be discussed 
below in more detail. 



IV. DISCUSSION 

A. Splitting versus Shifting of Bragg Peaks 

In this section we focus on the time region I (0 < 
t < Tsound) where we observe a splitting of the SRO 
Bragg peak for the thicker sample. We perform numer- 
ical lattice-dynamics calculations and use the results in 
order to simulate the transient XRD response of the SRO 
layer. 

It has recently been shown that a linear-chain model 
of masses and springs is well-suited to calculate the lat- 
tice dynamics triggered by the quasi-instantaneous ther- 
mal stress in SRO 32 . Moreover, the results of these cal- 
culations can be easily used to accurately simulate the 
transient x-ray response of such photoexcited nanolay- 
ered samples employing dynamical XRD theory^. We 
therefore apply this toolbox to the present case of laser- 
excited SRO layers of different thickness in order to elu- 
cidate the general features of the structural dynamics. 

The details of the linear-chain model simulations for 
the calculation of the photoinduced lattice dynamics of a 
thin SRO film on a STO substrate are given elsewhere^ 2 -. 
The essential ingredient is the assumption of an instan- 
taneous rise of an isotropic thermal stress in SRO due 
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FIG. 3. Calculated photo-induced strain field of the thicker 
SRO layer on STO for selected time delays. The laser excita- 
tion at t — launches coherent LA phonons which form sharp 
strain waves propagating through the layer and the substrate. 



to incoherently excited phonons at r = 0. In a recent 
UXRD study on hexagonal LuMnC>3 anisotropic elas- 
tic properties had to be accounted for which are due to 
the lower crystal symmetry compared to SRO^ 3 -. How- 
ever, the SRO unit cell only slightly deviates from cubic 
symmetry (pseudocubic) implying rather isotropic elas- 
tic properties^. The fact that the thermal stress in 
SRO builds up quasi-instantancously has been proven 
to be valid by several UXRD experiments on a few-ps 
timescale^Sr— . In fact, very tiny phase shifts of pho- 
toexcited coherent phonons modes in superlattices have 
been observed which evidence a crossover from a finite 
rise time (« 200 fs) to an instantaneous onset of the 
displacive thermal stress 2 ^. Nonetheless, as we discuss 
below, the assumption of an instantaneous driving force 
is sufficiently good for the dynamics considered in this re- 
port. In Section|TT]we mention that the very fast electron- 
phonon relaxation in SRO implies a thermal stress pro- 
file given by the exponential absorption of the pump 
light. Accordingly, we start our calculations by assuming 
an exponential thermal stress profile with a 1/e decay 
length defined by the optical penetration depth of SRO, 
£th = = 52 nm 2 !. 

The spatio-temporal strain field for the thicker SRO 
layer on STO resulting from the lattice dynamics cal- 
culations is shown in Fig. [3j The laser-induced thermal 
stress sets in at r = 0. The instantaneous rise of this driv- 
ing force coherently excites LA phonon modes up to very 
high frequencies. In total, the superposition of all co- 
herent phonons results in propagating strain wave fronts 
starting at the air-SRO and SRO-STO interfaces where 
the thermal stress is not balanced and exhibits large gra- 
dients. In case of the 94.8 nm SRO layer whose thickness 
is almost twice the optical penetration depth of 800 nm 
light, £op t , the dominating strain front is an expansion 



wave launched at the surface. The expansion wave start- 
ing at the SRO-STO interface has much lower amplitude 
but is still visible (blue line in Fig. EJ) . Thus the essen- 
tial feature inferred from the calculations for the thicker 
SRO layer is the generation of a propagating wave front 
dividing the SRO layer into an expanded sublayer near 
the surface whereas the remaining sublayer gets slightly 
compressed in total for < r < T S ound (green and red 
line in Fig. [3]). The thickness of these sublayers is grad- 
ually increasing and decreasing with time, respectively. 
The slight compression of the decreasing sublayer is ad- 
dressed in more detail below. At later times the coherent 
strain waves have propagated into the substrate forming 
a bipolar strain pulse^ 5 - leaving a stationarily expanded 
SRO layer. A very small acoustic mismatch of SRO and 
STO results in a negligible reflection coefficient of acous- 
tic waves at the SRO-STO interface 4 ^. That is, the strain 
waves do not travel back and forth several times inside 
the SRO layer which would lead to a breathing mode of 
this layer as observed in other material combinations or 
in free-standing film o 47 ' 48 . 

A closer look at the transient strain fields in Fig. [3] 
reveals that the propagating wavepackets are superim- 
posed by a peculiar fine structure which is most pro- 
nounced for the tensile component of the bipolar strain 
pulse. Since we solve the differential equations of the 
linear chain analytically these features are no numeri- 
cal artifacts^. In fact, these high-frequency oscillations 
are a characteristic feature of a discretized linear chain 
and do not occur in elastic continuum models* 5 -. They 
are a result of the fact that the motion on the linear 
chain is essentially initiated at the surface and interface 
due to the large gradients of the thermal stress at these 
points. In addition to the initial displacement of the out- 
ermost masses on the linear chain an oscillatory motion 
of these is launched 4 ^. A thorough description of this 
high-frequency component is out of the scope of this re- 
port and shall be given elsewhere. Note that these high- 
frequency modes require a very fast build-up of the ther- 
mal stress to be efficiently excited, i.e. a sufficiently short 
pump pulse and a fast electron-phonon relaxation. How- 
ever, in real SRO and STO crystals at room temperature 
the lifetime of acoustic phonons of such high frequency 
is expected to be very short due to anharmonic phonon- 
phonon scatterin g 31 ^ 50 ! 51 . 

In order to correlate the identified transient features 
of the photoinduccd lattice dynamics to the transient 
signals in UXRD experiments we calculate the diffrac- 
tion curves of the sample at each time step utilizing the 
transient lattice deformations presented in Fig. [3J Fig- 
ure H)(a) compares the numerical results with experimen- 
tal data similar to the data shown in Fig. [2(d) but at 
a higher pump fluence of 30 mJ/cm 2 and without sub- 
strate subtraction. The presented simulation accounts 
for the instrument function of the PXS and assumes a 
decay length of the thermal stress profile of £ t h = 44 nm 
to obtain a good match. This deviation from the optical 
penetration length gjt is addressed below. The Simula- 
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FIG. 4. (a) Comparison of the transient diffraction curves 
measured on the thicker SRO layer at the PXS (symbols) 
with simulations (solid lines). The experimental data were 
recorded at a pump fluence of 30 mj/cm 2 . Time delay in- 
creases from top to bottom and the curves are displaced for 
clarity. The simulations employ an exponential thermal stress 
profile with decay length f t h — 44 nm < £opt' (b) Transient 
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tion reproduces the UXRD data very precisely. In par- 
ticular, the splitting of the SRO Bragg peak and the in- 
dividual shifts as the peaks decay and grow are in perfect 
agreement. We can now correlate the features revealed 
by the UXRD data to the structural dynamics obtained 
from the numerical lattice-dynamics calculations (Fig. [3}. 
The initial rocking curve at r = represents the unex- 
cited SRO layer. As discussed above, the photoinduced 
thermal stress essentially launches a strain wave front 
at the surface which generates a gradually growing ex- 
panded sublayer. This rather thin expanded layer gives 
rise to the appearance of the broad and weak extra peak 
at lower angles. The peak grows as time delay increases 
since the thickness of this expanded sublayer increases. 
Simultaneously, the slightly negatively strained sublayer 



decreases in thickness leading to the gradual disappear- 
ance and broadening of the initial SRO peak. This ob- 
servation immediately verifies the very fast generation of 
the thermal stress which is required for the sharp division 
into differently strained sublayers^. 

The UXRD data also show a slight shifting of the indi- 
vidual peaks to higher angles as they rise and fall. As ev- 
idenced in Fig.rjja), the simulation also accurately repro- 
duces this behaviour. The reason of the transient shifts 
of both peaks can also be found in the calculated lattice 
dynamics and is related to the exponential profile of the 
thermal stress and of the resulting layer expansion. As 
can be inferred from Fig. [3J the expansion wave launched 
at the surface gradually imprints an exponential expan- 
sion profile into the SRO film. The average strain of 
the expanded sublayer thus decreases as time increases. 
Accordingly, the angular displacement of the related ris- 
ing Bragg peak is largest just after time zero and slowly 
decreases as time increases. The reason for the slight 
shift of the initial peak is a little more subtle. In simple 
words one may say that at early times the strongly ex- 
panding near-surface region of the SRO layer increasingly 
squeezes the less expanding near-interface region until 
the entire layer is expanded. In general, a gradient in the 
(instantaneous) thermal stress profile launches acoustic 
sound waves. In the considered case, the most domi- 
nant gradients appear at the surface and interface where 
the sharp strain waves are triggered. However, also the 
exponentially varying negative gradient inside the SRO 
film is responsible for a compressive strain wave compo- 
nent propagating towards the substrate. This compres- 
sive strain component is dominant over the tiny tensile 
strain wave from the SRO-STO interface (green and red 
line in Fig. [3]). Hence it causes an increasing compression 
of the corresponding sublayer as time increases which is 
represented by the slight shift of the initial SRO Bragg 
peak towards higher angles. At later times it turns out 
that the exponential gradient is also responsible for the 
corresponding shape of the bipolar strain pulse in the 
substrate (magenta and yellow line in Fig. [3J. 

The above explanations rely on the exponential depen- 
dence of the thermal stress and thus expansion profile 
which suggest that the particular shifting behaviour of 
the two peaks is sensitive to, e.g., the decay length £ t h 
of these profiles. To verify this we performed simula- 
tions with different values of this decay length. The final 
strain of the SRO layer and thus the Bragg angle of the 
SRO peak at late times is held constant. A comparison 
of the transient rocking curve at r = 6 ps with these 
simulations is shown in Fig. SJb). Indeed, the shorter 
the decay length the more shifted is the rising (decaying) 
peak due to a relatively larger average expansion (com- 
pression) of the near-surface (near-interface) sublayer at 
early times. Surprisingly, the best fit of the UXRD data 
shown in Fig. 0Ja) is achieved assuming a decay length 
of £th = 44 nm. As discussed in Section [Til the extraordi- 
nary fast electron-phonon relaxation in SRO should gen- 
erate a thermal stress and expansion profile having a de- 



7 



cay length £ t h equal to the optical absorption length £^pt . 
However, the decay length deduced from the UXRD data 
is significantly lower than the expected 52 nm. Even if 
one allowed for any ballistic or diffusive thermal trans- 
port before the energy is coupled into the lattice this 
result could not be explained since these processes tend 
to flatten out any gradients resulting in a longer decay 
length. 

There are two effects which could cause the unexpected 
steepness of the observed expansion profile. First, the 
optical penetration depth could possibly be different in 
the considered thin films as compared to the bulk ma- 
terial due to finite size effects and/or variations of the 
optical constants by stationary strains that are induced 
by epitax y 53 ' 54 . We performed spectroscopic ellipsome- 
try measurements on the thick SRO layer and found an 
optical penetration depth of — 48 nm at a wave- 
length of 800 nm. This value is indeed slightly smaller 
than the literature bulk value of 52 nm21 but still signifi- 
cantly larger than the lengthscale of the expansion profile 
deduced from the UXRD data. 

The second possible effect which could cause a steeper 
expansion profile is a weak temperature dependence of 
the Gruneisen parameter 7 describing the ratio of ther- 
mal expansion and deposited energy. In general, this 
parameter is nearly material-independent and shows al- 
most no temperature dependence, however, a slight in- 
crease of 7 with temperature can be observed in several 
materials^ 5 -. An increase of 7 with T would cause a larger 
expansion of the near-surface regions of the SRO layer rel- 
ative to the deposited energy density profile which would 
result in a steeper expansion profile as revealed by the 
UXRD data. 

In the following we briefly turn to the lattice dynamics 
of the thinner SRO layer whose thickness is much smaller 
than the optical penetration length of the 800 nm pump 
light. Here, we can extend the concept of differently 
strained sublayers. Since the deposited energy density is 
comparable near the surface and interface, respectively, 
the launched strain wave fronts are also similar in ampli- 
tude. Accordingly, one finds three sublayers of different 
strains which should in principle cause a more compli- 
cated splitting of the layer Bragg peak. However, due to 
the small thickness of the SRO layer the corresponding 
Bragg peak is very broad and thus prevents any splitting 
from being observed. The required strain amplitudes to 
observe a splitting for the very thin film would be too 
large. This complex splitting of Bragg peaks could pos- 
sibly be visible on thicker layers of materials having an 
accordingly larger optical penetration depth. 



B. Universal Features of Lattice Dynamics in SRO 
Thin Films 

In the previous section we discussed the qualitatively 
different UXRD signatures while the coherent strain 
waves pass through the layer once (0 < r < T SOU nd)- 



Region I Region II Region III 




Normalized Delay x / T 

sound 

FIG. 5. Transient average SRO strain for different ratios 



As already mentioned in section [TTT1 the transient shifts 
appear very similar for the two thin films with different 
thicknesses at times later than T soun( j (region II and III) . 
Moreover, the transient weighted shifts (average strain) 
are almost identical for all times. In the following we thus 
address the evolution of the average layer strain and how 
it depends on the ratio of the layer thickness d and the 
decay length of the thermal stress £th ■ 

In order to investigate the effect of a varying ratio 
d/^th, we calculate the spatio-temporal strain fields for 
the thick SRO layer (d fixed) at various values of £ t h- 
We then extract the transient average SRO strain nor- 
malized to the final strain at r > 2T sounc i and plot the 
results versus normalized time delay r/T soun( j in Fig. [5] 
We identify three regions of qualitatively different lattice 
dynamics which, however, show universal behaviour in- 
dependent of the ratio d/^th- In region I (0 < r < T soun d) 
the coherent strain waves cause the average SRO strain 
to monotonically increase to 150% of the final thermal 
expansion. This value of maximum strain does not de- 
pend on the ratio d/^th-, however, the precise evolution of 
the total SRO strain is influenced by the ratio. It changes 
from linear to jump-like as the ratio decreases. Similar 
observations hold for region II (T soun d < r < 2T soun( j) 
but here we find a monotonically decreasing strain until 
it reaches the final value determined by the pure thermal 
expansion due to the increased SRO temperature. In 
region III (r > 2T soun( j) the SRO strain remains quasi- 
constant until it relaxes back to zero via heat diffusion 
into the substrate on a nanosecond timescale (not in- 
cluded in the model and not shown). 

The quantitative difference of the total SRO strain evo- 
lution is simply explained by the different shapes of the 
strain waves launched at the surface and/or interface. In 
case of d/£ t h "C 1 (grey line in Fig. square-like strain 
pulses are launched at the surface and interface with com- 
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Region I Region II Region III 




Normalized Delay x / T 



FIG. 6. Transient normalized potential (solid lines) and ki- 
netic energy (dashed lines) of the coherent strain waves inside 
the SRO layer (red) and the STO substrate (black). The en- 
ergies were derived from the calculation presented in Fig. [3] 
Note that the potential and kinetic energy in the substrate 
are identical and thus the respective curves are overlayed. 
The black bullet indicates the point at which all normalized 
energies are identical. This point is independent of £ t h 

parable amplitudes which results in the piecewise linear 
evolution of the total strain. In contrast, if d/£th ^ 1 
(blue line in Fig. [5]) only a small portion underneath the 
surface of the SRO layer is excited which then rapidly ex- 
pands just after time zero. Moreover, a rather localized 
bipolar strain pulse is launched^ which has a vanishing 
integral strain as we discuss below. Around r = T soun( j 
the bipolar strain pulse leaves the SRO layer and causes 
the spike-like transient increase of the total SRO strain. 
Altogether, Figure E] evidences that the calculations pre- 
cisely predict the average SRO strain dynamics derived 
from the transient (weighted) shift of the SRO Bragg 
peaks shown in Fig. O In particular, the slight differ- 
ences of the red and blue data in Fig. [2 (SRO layers of 
different thickness) can thus mainly be attributed to the 
different ratios d/£th- 

The striking feature in Fig. [5] is that the normalized 
total strain in the SRO layer at r = T soun d and for 
t > 2T soun d (region III) is independent of the ratio d/£th- 
In particular, the maximum strain due to the coherent 
phonon dynamics is always 50% larger than the steady- 
state strain after 2T soun d due to incoherent phonons (i.e. 
heat). To give an explanation we employ the potential 
and kinetic energy of the coherent strain waves inside the 
SRO layer and the STO substrate, respectively, which are 
shown in Fig. El At r = the thermal stress is gener- 
ated quasi-instantaneously by the absorption of the laser 
pulse which increases the equilibrium distance of atoms 
in the lattice. Therefore, the thermal stress gives rise 
to an initial potential energy since the SRO layer is now 
compressed relative to the new equilibrium state. The 
kinetic energy is zero since no atomic motion has started 



at this point. As explained in section HV Al the thermal 
stress launches propagating strain waves inside the SRO 
layer which is evidenced by the increasing kinetic and de- 
creasing potential energy of SRO. The increase in both 
kinetic and potential energy inside the STO substrate is 
due to the corresponding strain pulse launched from the 
interface into the substrate (cf. Fig. . At r = T soun d all 
energies are identical and the key point is that this ob- 
servation is independent of the ratio G?/£th (not shown) 
which is indicated by the black bullet. In particular, the 
substrate-related energies reached half of their maximum 
values since the compressive part of the bipolar strain 
pulse has entered the substrate while the tensile part still 
remains in the SRO layer. This implies that the integral 
strain of the two individual parts of the bipolar strain 
pulse (tensile and compressive) is equal but has opposite 
signs. In contrast to the initially (r = 0) compressed 
state of the SRO layer relative to the new equilibrium 
state after 2T soun d, the coherent superposition of the lon- 
gitudinal strain waves at r = T S ound results in a tensile 
state relative to the new equilibrium (cf. Fig. At 
this point in time the potential energy of the SRO layer 
dropped to 1/4 of its initial value. Since the potential en- 
ergy generally is proportional to the square of the strain 
we conclude that the relative expansion at r = T soun d is 
1/2 of the initial relative compression. Hence the max- 
imum absolute strain relative to the unexcited state at 
t = T soun d has to be 150% of the equilibrium (thermal) 
strain after 2T soun d. 



As mentioned in section IIV Al the perfect acoustic 
matching of SRO and STO prevents the coherent strain 
waves from being (partially) reflected back into the SRO 
layer. If there is a significant acoustic mismatch the 
strain waves would travel back and forth until all vibra- 
tional energy has been transferred to the substrate. This 
would result in an oscillatory behaviour of the total SRO 
strain as was observed by UXRD in a photoexcited Gold 
layer on a Mica substrate^. The limiting case of such an 
acoustic breathing of a metal film is represented by free- 
standing films where the amplitude decay of the acoustic 
waves is merely given by internal damping effects. Such 
breathing of free-standing Al films has been studied by 
femtosecond electron diffraction^. Ideally, the maximum 
strain at t = Tsound due to the coherent structural dy- 
namics in an acoustically decoupled photoexcited thin 
film is 100% above the steady-state thermal expansion 
which is related to the displacive nature of the excitation 
mechanismiS,!^. Any type of energy loss — be it coherent 
by transmission losses into the substrate or incoherent 
by internal damping and/or scattering — leads to a de- 
crease of the maximum coherent strain at r = T S ound- In 
the present case of a semi-infinite substrate with perfect 
acoustic matching the coherent strain maximum of the 
metal layer is reduced by a factor of 2 since a part of the 
initial potential energy is transferred to the substrate in 
form of the compressive half of a bipolar strain pulse. 
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V. CONCLUSION 

In this report we address the issue of the one- 
dimensional coherent structural dynamics in thin films of 
metallic SrRu0 3 (SRO) on a dielectric substrate SrTi0 3 
which are triggered by the absorption of ultrafast opti- 
cal laser pulses. We show experimental results of UXRD 
experiments which probe the photoexcited coherent lat- 
tice dynamics of two SRO layers with different thick- 
ness. The observed changes of the respective SRO Bragg 
peaks exhibit qualitatively different features. In particu- 
lar, the Bragg peak of the thicker SRO layer shows a tran- 
sient splitting into two separated reflections which evi- 
dences the propagation of sharp longitudinal strain waves 
through the metal layer. These coherent wavepackets are 
caused by the exceptionally fast electron-phonon relax- 
ation previously identified in SRO^ 4 . Numerical models 
accounting for the photoinduced structural dynamics and 
dynamical XRD show excellent agreement assuming a 
surprisingly small decay length of the thermal stress pro- 
file of 44 nm. This deviation from the measured optical 



penetration depth in these thin films may be attributed 
to a temperature-dependent Gruneisen parameter. Fi- 
nally, we analyze the UXRD-calibrated coherent phonon 
dynamics in detail using the numerical simulations. The 
UXRD features can unambiguously be related to the pre- 
cise structural dynamics. We discuss the features specific 
to the layers of different thickness (Bragg peak splitting 
versus continuous shift) and identify a universal evolu- 
tion of the total strain of a photoexcited thin metal film 
of arbitrary thickness on a semi-infinite substrate. Our 
work gives a very precise and UXRD calibrated picture 
of the acoustic deformation of laser-heated metal layers 
on a supporting substrate and carefully relates the tran- 
sient structural features to UXRD signatures. We believe 
that this work is very valuable for the quantitative inter- 
pretation of time-resolved scattering experiments on the 
complex photoinduced structural dynamics of crystalline 
nanolayered structures such as thin films, multilayers and 
superlattices. 

We thank the BMBF for funding the project via grant 
No. 05K10IP1 and the DFG via grant No. BA2281/3-1. 
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